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The standard model extended by three vector-like families with masses of order 1 TeV - 100 TeV 
allows for unification of gauge couplings. The values of gauge couplings at the electroweak scale are 
highly insensitive to fundamental parameters. The grand unification scale is large enough to avoid 
the problem with fast proton decay. The electroweak minimum of the Higgs potential is stable. 



Introduction. Our current understanding of elemen- 
tary particles and their interactions is described by the 
standard model of particle physics (SM). It contains 
many free parameters, namely masses and mixing of 
quarks and leptons of the three families, the Higgs boson 
mass, and most importantly, the values of three gauge 
couplings ai,2,3 that determine the strengths of electro- 
magnetic, weak, and strong interactions. One of the main 
goals in particle physics is to understand values of these 
free parameters from basic principles. 

Among the most elegant approaches is the idea of 
grand unification in which three gauge couplings q;i_2,3, 
corresponding to three different symmetries of the SM, 
originate from a single gauge coupling associated with the 
symmetry of a grand unified theory (GUT) 1 . However, 
in the standard model the gauge couplings do not unify. 
Although through renormalization group (RG) evolution 
they run to comparable values at about 10^* GeV the 
mismatch is too large to be accounted for by GUT scale 
threshold corrections. In addition, the value of the uni- 
fication scale is too low to satisfy limits on proton de- 
cay [1]. Thus, in order to realize this idea, the SM has to 
be extended by additional particles. 

In this letter we show that extending the standard 
model by three complete vector-like families (3VF) with 
masses of order 1 TeV - 100 TeV allows for unification of 
gauge couplings. Predictions for gauge couplings at the 
electro-weak (EW) scale are highly insensitive to funda- 
mental parameters, and ratios of observed values are to 
a large extent understood from the particle spectrum it- 
self. The GUT scale can be sufficiently large to avoid the 
problem with fast proton decay, thus resurrecting simple 
non-supersymmetric GUT models. 

The way this scenario works can be summarized in few 
steps. First, extra 3VF make all gauge couplings asymp- 
totically divergent which opens a possibility for a unifi- 
cation with large (but still perturbative) unified gauge 
coupling. Consequently, in RG evolution to lower ener- 
gies gauge couplings run to the infrared fixed point. Sec- 
ond, the ratios of gauge couplings close to the infrared 
fixed point depend mostly on the particle content of the 
theory and they happen to be not far from the observed 
values. Finally, the discrepancies between values of gauge 
couplings predicted from closeness to infrared fixed point 
and corresponding observed values can be fully explained 



by threshold effects from masses of particles originating 
from 3VF. We note, that the first part is very similar to 
attempts to explain observed values of gauge couplings 
from infrared fixed point with 8 to 10 chiral families [3] 
(see also Refs. [HIS]) before the number of chiral families 
and values of gauge couplings were tightly constrained. 

Allowing arbitrary new particles provides many possi- 
bilities for gauge coupling unification. However, with ar- 
bitrary new particles the predictive power is typically lost 
and a GUT embedding of such models is more compli- 
cated. On the other hand, the addition of complete fam- 
ilies, either chiral (only left-handed or only right-handed 
particles with given quantum numbers) or vector-like (for 
each left-handed particle there is a right-handed particle 
with the same quantum numbers) represents some of the 
simplest extensions of the SM (especially since we have no 
understanding of why there should be just three families 
in nature), that can be easily embedded into GUTs. We 
will also argue that the predictive power of this scenario 
is comparable to minimal supersymmetric unification. 

While the existence of extra chiral families is unlikely 
given the current experimental situation, the vector-like 
families are poorly constrained ^6j. This is mostly be- 
cause a vector-like pair of fermions can have a mass 
which is not related to their couplings to the Higgs bo- 
son. Depending on the dominant decay mode, the lim- 
its on new vector-like fermions range from ^ 100 GcV 
to ^ 500 GeV. Consequently, there is a vast literature 
on models that contain extra vector-like fermions near 
the TeV scale. Examples include attempts to explain 
the anomaly in the forward-backward asymmetry of the 
b-quark |S1 [TO] , the muon g-2 anomaly [8^ , and studies 
of effects of extra VFs in supersymmetric theories on 
gauge coupling unification |11L 112) . spectrum of super- 
partners 13 , or the Higgs mass [7], among many others. 
Adding three vector-like families. Consider the one- 
loop RG equations (RGEs) for three gauge couplings : 

(1) 

where t — In Q/Qo with Q representing the energy scale 
at which gauge couplings are evaluated, and the beta 
function coefficients in the SM are given by 
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where is the number of famihes. For n/ = 3 we get the 
usual resuh, bi — (41/10, —19/6, —7), while with extra 
3VF, n/ = 3 + 2x3 = 9(a vector-like partner contributes 
in the same way), we find 6, = (121/10, 29/6, +1). Thus 
in the SM-I-3VF all three couplings are asymptotically 
divergent. 

An example of the RG evolution of gauge couplings 
in the SM-I-3VF with unified gauge coupling ac = 0.3 
at Mg = 2 X 10^^ GeV is given in Fig. [ij The Mq is 
chosen such that the limits on proton decay are satis- 
fied. Our choice also coincides with the scale of min- 
imal supersymmetric unification and thus allows for a 
direct comparison. The evolution of gauge couplings in 
the SM, that will be useful for discussion of threshold 
corrections, is showed by dashed lines. As input we use 
the experimental central values of a~^\j{Mz) = 127.916 
and sin^^vK = 0.2313, which are related to ai,2(-^^z) 
through sin^ 0w = «'/ (a2 +a') and asM = 0-2 sin^ dwi 
where, using the SU{5) normalization of hypercharge, 
a' = (3/5)ai; then a3{Mz) = 0.1184; and the top quark 
mass mt — 173.2 GeV. These together with current ex- 
perimental uncertainties can be found in Ref. [5]. We 
neglect Yukawa couplings of other fermions in the SM, 
and for simplicity, we also assume that Yukawa couplings 
of 3VF are negligible. We further assume the Higgs bo- 
son mass m/i = 125 GeV which is in the middle of the 
currently allowed range [13]. In all numerical results we 
use full two loop RGEs [15]. All particles with masses 
above Mz are integrated out at their mass scale. We in- 
clude one-loop matching corrections for and ruh |16| . 
However, in order to understand results, approximate an- 
alytic formulas will be sufficient. 

The one-loop RGEs, which are very good approxima- 
tions for ai_2, can be easily solved, and gauge couplings 
at the EW scale can be written in terms of the GUT scale, 
Mc, and the unified gauge coupling, (^^(A/g) ~ ac- 
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where Ti are the threshold corrections that depend on 
masses of the extra vector-like fermions. They can be 
approximated by the leading logarithmic corrections: 
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where bj is the contribution of a given fermion /, 
with mass Mf, to the corresponding beta function co- 
efficient [in]. Neglecting these corrections for a mo- 
ment, we find that for large enough unification scale 
and sizable but still perturbative unified gauge cou- 
pling the first term in Eq. ([s]) dominates, a~^{Mz) — 
{bi/2n)\n{MQ/Mz) , and the ratios of gauge couplings 
are completely fixed by ratios of beta function coeffi- 
cients, ai{Mz) /aj{Mz) — bj/bi. Alternatively, in this 0- 
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FIG. 1: RG running of gauge couplings: as (top solid line), 
a2 (middle solid line), and ai (bottom solid line) in the SM 
extended by three vector-like families for ac ~ 0.3 at Mq = 
2 X 10^^ GeV. Dashed lines in the same order show running 
of gauge couplings in the SM. Masses of 3VF are neglected in 
the top plot, and fixed to 10 TeV in the bottom plot. 



th order approximation, we obtain a parameter-free pre- 
diction for sin^ 9w'- 
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where b' = (5/3) 61, and the numerical value corresponds 
to the SM-I-3VF. It is identical to the one obtained as- 
suming 9 chiral families [3H5]. Although it is not a per- 
fect match to the measured value, we will see that the 
discrepancy can be easily accommodated by taking into 
account threshold corrections, T^, and finite value of ac- 
For aa, the one-loop RGE, Eq. ([T]), is not a good 
approximation because of the accidentally small 63 fac- 
tor. The two-loop contribution to the beta function, 
ai^s/STT^, where B3 = -102 + (76/3)n/ = 126 for 
SM-I-3VF [TS], is larger than the one- loop contribution 
for 0:3 > 0.1 (the coupling is perturbative, the 3-loop 
contribution is only a small correction for a3 ~ 0.1). 
Neglecting the one-loop contribution and threshold cor- 
rections, the RGE for 03 can be easily solved: 
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Again, for our Mq and og, the log term dominates. 
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Numerically, for aEM{Mz) = 1/127.916, it predicts 
a3{Mz) = 0.099. The 1-loop contribution can be added 
as an expansion in e = Anbs/Bs: a^^Mz) — >■ a3{Mz)/K, 
where K=l + (l/3)(e/a3(Mz)) - (l/12)(e/a3(Mz))2 + 
.... It represents ~ 26% correction and leads to 
a3{Mz) — 0.073 which is close to 0.066 obtained using 
full two-loop RGEs, see Fig. [l] (top). The discrepancy 
from the observed value can be accommodated by the 
threshold correction, T3, which is well approximated by 
Eq. (j4]) and can be added to a^'^{Mz) as in Eq. i3l. 
Threshold corrections. Assuming a common mass, 
MvF, of all particles from extra 3VF, the threshold cor- 
rections, given in Eq. Q, have the same value in the lead- 
ing log approximation, Ti = {A/Tr)\n{MvF/Mz), which 
is ^ 6 for MvF = 10 TeV. These modify gauge couplings 
at the EW scale by factor 1 -I- ai{Mz)Ti, which repre- 
sents ~ 10%, ~ 20%, and ~ 40% corrections to gauge 
couplings ai, a2, and respectively. For our example 
with ac = 0.3 at Mq = 2 x 10^*^ GeV, the common 
mass MvF — 10 TeV leads to gauge couplings that are 
in agreement with experimental values with better than 
8% precision. The effect of Mvf = 10 TeV is shown 
in Fig. [1] (bottom). This prediction is highly insensitive 
to the value of ac- From Eqs. ^ and ^ we can see 
that any aa ^ 0.3 contributes less than ~ 10% to gauge 
couplings at the EW scale. Its exact value is important 
only for precise prediction for gauge couplings as are ex- 
act masses of all particles from extra 3VF in case they 
are split. Thus the only two relevant parameters are Mq 
and Mvf, and with better than 8% agreement of pre- 
dicted values of three gauge couplings with experimental 
values, the predictive power of this scenario is compara- 
ble to that of minimal supersymmetric grand unification 
(which predicts ^ 10% too large a3{Mz); this discrep- 
ancy is usually explained by 3% - 4% GUT scale threshold 
corrections, or by splitting superpartner masses [1]). 
Sensitivity to fundamental parameters. Since ac ^ 
0.3 contributes less than ^ 10% to the EW scale values 
of gauge couplings, changing it by a factor of 2 does not 
modify predicted values of gauge couplings by more than 
^ 10%. This can be contrasted with high sensitivity in 
the minimal supersymmetric unification in which a^^ ~ 
24, and ~ 3 times larger value of a3{Mz) compared to 
aa is the result of a cancellation between the \n{MG /Mz) 
and a^^ terms in Eq. Thus a given variation, x%, 
of aa results in ~ 3a;% variation of a3{Mz). 

Sensitivity of predicted gauge couplings to Mq and 
Mvf is also very small. Changing Mq by a factor of 
two changes \n{MG/Mz) term in Eqs. ^ and ^ and 
thus the weak scale values of gauge couplings by ^ 2%. 
In minimal supersymmetric unification the same change 
would result in ~ 6% change in a3{Mz) again as a result 
of the above mentioned cancellation. 

Finally, changing Mvf by a factor of two changes 
threshold corrections by ~ 15% which in turn changes 
EW scale values of gauge couplings by no more than 



6% (threshold corrections represent ^ 40% of a3{Mz))- 
Thus, changing any of the fundamental parameters by a 
factor of 2 does not modify predicted values of gauge cou- 
plings by more than ^10%. 

Realistic example. In order to obtain gauge couplings 
within experimental uncertainties, the masses of parti- 
cles from extra 3VF must be split. This is indicated in 
Fig. [1] (top) by close but not identical scales at which 
the RG evolutions of gauge couplings in the SM and 
SM-I-3VF cross. These crossing scales are Mi ~ 100 
TeV, M2 ~ 1 TeV, and M3 ~ 10 TeV for ai, az, and a^, 
respectively. Thus any spectrum of particles from 3VF 
for which the threshold corrections given in Eq. Q equal 
to Ti = {4:/TT)ln{Mi/Mz), as if all particles charged un- 
der given gauge symmetry had common mass Mi, will 
reproduce the measured values of gauge couplings. 

There are many solutions available. For simplicity, we 
present one that assumes universal masses for all vector- 
like pairs with the same quantum numbers. We use the 
same names for vector-like paris as for quarks and lep- 
tons in the SM, e.g. mq is the common mass of vector- 
like quark doublets. For ac = 0.3 at Mq = 2 x 10^^ 
GeV a fully realistic example is obtained with mq = 500 
GeV, TOL = 95 TeV, mu = 220 TeV, = 180 TeV, 
ruE = 250 TeV, and is shown in Fig. [2] (top). These pa- 
rameters, specified to two significant figures, reproduce 
measured values of gauge couplings with much better 
than 0.1% precision and thus within small fractions of 
current experimental uncertainties. This is another way 
to demonstrate high insensitivity of gauge couplings to 
fundamental parameters. 

For completeness, evolutions of the top Yukawa and 
the Higgs quartic couplings are given in Fig. [2] (bottom). 
The Higgs quartic coupling remains positive all the way 
to the GUT scale and thus the electroweak minimum of 
the Higgs potential is stable. This result holds for any 
Higgs boson mass in the currently allowed range [T3]. 
Discussion and conclusions. Since the effect of inte- 
grating out vector-like families can only increase gauge 
couplings at the EW scale, see Eqs. pp] ) and Fig. [T] 
values of predicted couplings without considering mass 
effect of 3VF have to be smaller than the measured val- 
ues. In other words, the lines representing evolution of 
gauge couplings in the SM-I-3VF have to cross those of 
the SM before they reach the Mz scale. This leads to a 
lower limit on the GUT scale. For aa = 0.3 the lower 
limit is Mg 1 x 10^^ GeV. There is no upper limit on 
Mg, and it can be identified with the string scale or the 
planck scale. Increasing Mq however increases the aver- 
age mass scale of 3VF, that can be inferred from Eq. ([3| 
or Fig. [1] and also increases the sensitivity of predicted 
gauge couplings to Mvf- We checked numerically that 
increasing AIg also requires larger splitting of masses of 
particles from 3VF. 

Adding more vector-like families leads to faster run- 
ning, and correspondingly the masses of VF must be 
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FIG. 2: Top: the same as in Fig. [T] with masses of three pairs 
of vector-hke fermions fixed to: mq = 500 GeV, ttil = 95 
TeV, mu = 220 TeV, mo = 180 TeV, and ruE = 250 TeV. 
Bottom: renormalization group running of the top Yukawa 
coupling, At, and the Higgs quartic coupling, A^. 



Yukawa couplings to a singlet scalar field, or the extra 
3VF might be charged under additional symmetries. 

Similar unification mechanism is expected to work in 
other models that do not modify RGBs for gauge cou- 
plings dramatically, for example the two Higgs doublet 
model or other simple extensions of the SM. 

Another possible interpretation that this scenario of- 
fers is a non-perturbative unification at a high scale, as 
envisioned in Ref. 3J. However, it also allows an in- 
terpretation without any GUTs: gauge couplings might 
originate from random large values at a high scale. This 
follows from the fact that ai{Mz) very weakly depend on 
ai{MG), and thus whether gauge couplings really unify 
or not is not important for understanding the EW scale 
values of gauge couplings. However, this would not pro- 
vide an explanation of quantum numbers of SM particles 
which is perhaps the main virtue of simple GUTs. 

The scenario we present does not offer an explanation 
of the hierarchy between the EW scale and the GUT 
scale. This might not be a disadvantage, since we have 
not discovered any signs of new physics that would solve 
the hierarchy problem yet. 
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larger or the GUT scale must be significantly lowered. 
Lowering the GUT scale is disfavored by proton decay 
limits, and increasing masses of VF (for example, with 
4VF the average mass scale of VF increases to ^ 10^ GeV 
for Mq = 2 X 10^® GeV) leads to larger sensitivity of 
predicted gauge couplings to My p. Thus, the SM-f-3VF 
represents the minimal and the most attractive scenario 
for insensitive unification. 

Although there are many possible arrangements for 
masses of particles from 3VF that lead to precise pre- 
dictions for gauge couplings, the existence of a solution 
with relatively small splitting between masses and their 
proximity to the EW scale is not trivial. Requiring the 
smallest sensitivity to masses of particles from 3VF sug- 
gests that some of the extra fermions might be within 
the reach of the LHC. Since a2 has the lowest cross- 
ing scale, this applies especially to quark doublets Ql,r 
and lepton doublets Li^jj. However, we checked numer- 
ically that breaking mass universality among particles 
with the same quantum numbers, each particle separately 
can be at current experimental limits. A lot of freedom 
in masses of extra VF indicates that it should be possible 
to combine explanations of various anomalies mentioned 
in the introduction with insensitive unification. It might 
also allow for simple GUT scale boundary conditions for 
these masses (this freedom is further enhanced by possi- 
ble Yukawa couplings that change their RG evolution). 
Alternatively, the masses of 3VF might originate from 
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For a review and references, see the section on grand 
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